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The NASA Fast Scattering Code (FSC) is a versatile noise prediction program designed 
to conduct aeroacoustic noise reduction studies. The equivalent source method is used to 
solve an exterior Helmholtz boundary value problem with an impedance type boundary 
condition. The solution process in FSC v2.0 requires direct manipulation of a large, dense 
system of linear equations, limiting the applicability of the code to small scales and/or 
moderate excitation frequencies. Recent advances in the use of Fast Multipole Methods 
(FMM) for solving scattering problems, coupled with sparse linear algebra techniques, 
suggest that a substantial reduction in computer resource utilization over conventional 
solution approaches can be obtained. Implementation of the single level FMM (SLFMM) 
and a variant of the Conjugate Gradient Method (CGM) into the FSC is discussed in this 
paper. The culmination of this effort, FSC v3.0, was used to generate solutions for three 
configurations of interest. Benchmarking against previously obtained simulations indicate 
that a twenty-fold reduction in computational memory and up to a four-fold reduction in 
computer time have been achieved on a single processor. 
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Nomenclature 

FSC complex scattering matrix, A e c NxM 
Adjoint matrix, A* e C Mx7V 

The linear space of complex matrices with N rows and M columns 
The linear space of complex column vectors of length N 
Speed of sound 

Spherical Bessel function of order V 

Second kind spherical Hankel function of order V 

Excitation wave number 

Number of receiver boxes for SLFMM discretization 
Number of source boxes for SLFMM discretization 
Unit surface normal vector directed into the fluid 
Number of collocation points on scattering surface(s) 

Number of equivalent sources on source surface(s) 

Incident acoustic pressure 
Scattered acoustic pressure 

associated Legendre function of degree v and order jU 
Spherical harmonic expansion truncation number 

Receiver boxes 

Regular basis functions for harmonic expansions 
Scattering geometry surface area 
Source boxes 
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fcE 

= 

Source box enveloping sphere interiors 

kt 

= 

Source box enveloping sphere exteriors 

s; 

= 

Singular basis functions for harmonic expansions 

— / 

U ino 

= 

Incident acoustic velocity 

U scat 

= 

Scattered acoustic velocity 

(x,t) 

= 

Space-time coordinates 

Mi, 

= 

Collocation points on scattering surfaces 

{*• f 

L m J/w=l 

= 

Receiver box centers 


= 

Source points on source surfaces 

fcE 

= 

Source box centers 

y m 

± v 

= 

Spherical harmonics 

a 

= 

Surface normal admittance 

ip,o,<p) 

= 

Spherical coordinates 

Po 

= 

Fluid density 

CD 

= 

Excitation frequency 


I. Introduction 

Successful aircraft deployment requires an accurate assessment of its impact on community noise. Preceding any 
experimental study of noise abatement techniques, engineers rely on physically accurate computational methods to 
reduce the set of parameters to be investigated. The NASA Fast Scattering Code (FSC), developed at the Langley 
Research Center (LaRC), is a versatile noise prediction program for evaluating noise reduction concepts such as the 
shielding effects of scattering geometry and engine nacelle placement, or sound absorption through liner 
optimization. The code calculates the time harmonic acoustic field produced by the interaction of known incident 
sound with aero-surfaces, of arbitrary shape and surface impedance, immersed in a potential flow. Researchers have 
applied the FSC to a variety of aeroacoustic simulations involving fixed-wing and rotary-wing configurations, as 
detailed in references 1-5. 


The FSC uses the equivalent source method (ESM) to solve an exterior three-dimensional Helmholtz boundary 
value problem (BVP) by expanding the scattered acoustic pressure field into a series of point monopoles or dipoles. 
The point sources are located on a fictitious surface located interior to the actual scattering surface. Users input 
surface geometry, background fluid properties, incident fields for both acoustic pressure and velocity, and operating 
conditions. The code then calculates the scattered acoustic pressure field at user specified locations in space and 
time. Details on the solution methodology and program usage are given in references 6 and 7. 


The use of boundary techniques - such as the ESM, finite element methods, boundary integral methods, or 
boundary element methods - for solving BVPs involves discretization of the scattering surfaces and the solution of a 
large (usually) dense system of linear equations. For the ESM, an over-determined linear system is produced where 
the number of unknown equivalent monopole source strengths is approximately equal to the product of scattering 
surface area ( S) and the square of the incident sound excitation wave number ( k)\ 


N s ~Sk 2 


( 1 ) 


As implemented in the FSC, the ESM requires approximately one-third the memory of the other BVP solution 
methods and about one-ninth the computational time (accounting for the adjective Fast in FSC). Even with these 
computational gains, FSC version 2.0 predictions using current single-processor workstation technology with less 
than 30 GB of shared memory are limited to conditions where the number of equivalent sources (N s ) does not 
exceed about 25000. Simulations for high frequency/large scattering surface combinations, such as a full scale 
commercial transport, yield large linear systems with considerably more unknowns that can only be solved with 
advanced computer systems or by employing numerical algorithms that drastically reduce computational resource 
utilization. 
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Recent advances in the use of Fast Multipole Methods (FMM) for solving scattering problems coupled with 
sparse linear algebra techniques, suggest a substantial reduction in the consumption of computer resources relative 
to dense matrix based boundary methods 8 . Implementation of FMM into the FSC is described in the present paper. 
In the resulting computer program, FSC v3.0, the dense, linear least squares solution approach of FSC v2.0 is 
replaced with an iterative conjugate gradient method (CGM) that features FMM based fast field calculations and fast 
matrix- vector multiplication algorithms. Numerical scattering studies are presented that show an order-of-magnitude 
reduction in computer memory and computational time. This results in the ability to quadruple the number of 
unknowns in the linear system to be solved, which equates to doubling the maximum frequency achievable with 
FSC v2.0 for a given problem. Discussions in the remainder of the paper are grouped into two categories: 1) FSC 
computational advances, and 2) aeroacoustic applications including high frequency numerical simulations for full 
scale nacelle and commercial transport configurations. 


II. FSC Computational Advances 

Technical aspects of fast multipole methods as they pertain to the FSC version 3.0 are discussed in this section. 
First, the Single-Level Fast Multipole Method (SLFMM) is introduced for the rapid calculation of the scattered field 
produced by a collection of monopoles. Next, memory reduction algorithms are developed for solving the FSC 
linear system using conjugate gradient methods in which the SLFMM is used for fast vector-matrix multiplications 
at each iteration of the CGM. Much of the terminology and notation used for the FMM analysis presented here was 
taken from reference 8. This comprehensive text contains detailed analyses of the FMM concepts and a wealth of 
references that follow the historical developments of the subject and its research applications. 


A. Single Level Fast Multipole Method 


For the monopole variant of the ESM, the scattered acoustic pressure at any spatial location x is written as a 
point monopole superposition (e + ^time dependence is assumed): 


/C,(*)=Z 


N s e -ik\x-y m \ 


47r\. 


x-y. 


( 2 ) 


In Eq. (2), {y m } N m s =l are the spatial locations of the source points and the complex numbers {c m are the source 
strengths. The source points are distributed inside the scattering geometry and the expansion strengths are tuned to 
satisfy a general impedance type boundary condition at the N c distinct collocation points {x n that lie on the 
scattering surfaces. For FSC applications N c ~ 3N S . The resulting over-determined linear system is solved using 
least squares methods as discussed in the next sub-section. 

As it pertains to the ESM, the basic concept of the FMM is to replace the field due to a single point monopole 
located at y by an infinite series of spherical harmonics in which the source point and field (or receiver) point are 
separated. The fundamental identity of the process is given by Eq. (3) and can be derived by applying the method of 
separation of variables to the three-dimensional Helmholtz differential equation in spherical coordinates. 


- **Z Z s v m ( y)K (*) |*| < 1*1 

v=0 ju--v 

< 

- **Z Z R y M (y) s v (*) 1*1 < 1*1 

v=0 jU=-v 


( 3 ) 


In (3), 

K{l )= j v {kp)Y?(e,(p) (4a) 

S^)=h[ T) (kp)Y l u (e,(p) (4b) 

are referred to as regular and singular basis functions, respectively, where (p,0,(p) are spherical coordinates 
corresponding to the Cartesian coordinates ^ , j v is the first kind spherical Bessel function of order v , h ( v 2) is the 
second kind spherical Hankel function of order v , and the spherical harmonics are given by the formulas 
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( 5 ) 


with denoting the associated Legendre function of degree v and order \jLl\ . The singular and regular 

terminology refers to the behavior of the spherical Bessel functions at p = 0 . Fast multipole algorithms feature 
truncated spherical harmonic expansions of both the singular and regular type for the approximation of the scattered 
acoustic field of equation (2). Direct evaluation of Eq. (2) at N receiver points requires 0(NN s ) mathematical 
operations. By achieving a large degree of source-receiver point separation, the FMM can significantly reduce the 
computational complexity of field calculations as we now demonstrate. 

Suppose all source points lie inside some sphere of finite radius. Let B + and B~ denote the exterior and interior 
regions of the sphere, respectively; then combining Eqs. (2) and (3) and formally interchanging the order of 
summation yields the so-called middleman algorithm (MA) for calculating the scattered field: 


PscaXx) = 


ZZW(x) xeB~ 

v=0 ju=—v 
v=0 ju=—v 


( 6 ) 


Where 


C: =-ikf j c m S-/(y m ) (7a) 

m = 0 

D^-ik^R/iyJ (7b) 

m = 0 


If for some prescribed truncation error tolerance £ , the expansions in Eq. (6) converge after the first v = q - 1 
terms, then the number of mathematical operations required to evaluate Eq. (6) at N receiver points is o(q 2 N) and 
o(q 2 N s ) operations to compute Eq. (7). Furthermore, if q 4 2 «N S , then the total mathematical complexity is 
0(N + N s ) compared to 0(NN s ) for direct evaluation using Eq. (2). The computational gains are entirely 
attributable to the separation of source and receiver interaction. 

For general scattering problems, the middleman algorithm fails because the sphere that envelops all sources is 
too large for the rapid convergence of Eq. (6). The SLFMM algorithm applies the MA to smaller spatial regions and 
has been implemented into the FSC. The process begins with the discretization of the computational domain into 
cubical boxes and the partitioning of source and receiver points into the boxes. The box size depends on the 
excitation frequency and size of the computational domain. Note that some boxes may contain both source and 
receiver points. Figure 1 shows the partitioning scheme for a futuristic large-scale commercial transport with three 
engine nacelles (SAX-40) 9 . Point monopoles are distributed over the entire source surface which lies slightly interior 
to the actual scattering surface. For this demonstration, it is desired to compute the scattered acoustic field on a plane 
of discrete observer locations bisecting the airframe as shown in the figure. 
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Receiver box 
partitioning 


Figure 1. Source and receiver point partitioning for the SAX-40 configuration, f = 400 Hz. 


Let {S n denote the collection of non-empty source boxes with box centers feE . and RXi the non- 
empty receiver boxes with centers We use the notation and S~ for the exterior and interior, 

respectively, of the smallest sphere that totally contains the n-th source box. A similar notation is used for the 
receiver box array. 

/T 

The radius of the sphere is given by d = d cube , where d cube is the length of the side of a cube. The term well- 

separated is used to describe two boxes whose centers are greater than 2d apart. The superposition given by Eq. (2) 
is distributed across the K s source boxes 


K s K s -ik\x-y a \ 

pL (*) = i>» (*) = £ z | 

n = 1 n = 1 y a eS n \ x Joe 


( 8 ) 


and evaluated according to the SLFMM algorithm as follows: 

Step 1) For n = , compute singular expansion coefficients for source box n: 

D? = -ik (?« - fn ) ( 9 ) 

y a ^ s n 

Step 2) For each receiver box R m , m = \,...,K R calculate t nm =x* m ~y* n and determine if source and receiver boxes 
are well-separated. If |F wiw | > 2d , then the boxes are well-separated and the regular expansion is used to compute the 
field for each point in the receiver box: 

/>;(*)=! I wM:) (io) 

v=0 ju=-v 

Where 

C:=-ik£c a S/{y a -? m ) ( 11 ) 

ya G ^n 
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If y nm \ < 2d , then source and receiver boxes are not well- separated and the singular expansion is used to calculate 
the field for receiver points outside the source box sphere, or the direct summation is used for receiver points that lie 
inside the source box sphere: 


Pnd) 



e ~ ik\x-y a \ 

*-y a 


xsS; n R m 


tes;nR m 


( 12 ) 


For large-scale scattering problems, most source-receiver box combinations are well-separated and the 
computational cost is dominated by the evaluation of Eqs. (10) and (11). Direct evaluation of the monopole sum in 
Eq. (12) occurs infrequently, and for some applications considered here accounts for less than 5% of the 
computational time. If the number of receiver points is of the same order of magnitude as source points, then 
analysis shows 8 that the computational complexity of the SLFMM is olq^Np ] compared to o(tv 2 ) for 100% 
direct evaluation. For q 2 «N S the savings are considerable. Further speed-up is possible by employing fast 
translation methods for the evaluation of the regular expansion coefficients of Eq. (11) given the singular expansion 
coefficients of Eq. (9). As discussed in detail in chapters 3, 7, and 8 of reference 8, fast translation algorithms reduce 
the SLFMM complexity to O^qN^^J • 

The SLFMM algorithm presented in Eqs. (8) through (12) has been implemented into the FSC v3.0 for 
predicting the scattered acoustic pressure field and its gradient and applied to known solutions from FSC v2.0 for 
comparison as reported in section III. Pressure gradient calculations require application of the gradient operator to 
the singular and regular basis functions of Eq. (4). Analytical details of this process are lengthy and the resulting 
recursion relationships can be found in chapter 2 (section 2) of reference 8. The SLFMM has also been deployed in 
the FSC v3.0 for fast matrix- vector multiplication during the linear algebra phase of the FSC process, as described 
below. 

B. Linear Algebra Considerations 

The FSC least squares matrix A e c N ° xNs is obtained by requiring the expansion given in Eq. (2) to satisfy an 
impedance type boundary condition at the N c collocation points {x}^ c =l on the scattering surface(s): 


AyW PoCok 1 


■^g im +(X,g, m j = \...,N c m = \,...,N s N C >N S 


( 1 ) 


Where 


Sim = 


4nR. 


R im = Xf-y m 


(14) 


In Eq. (13), hj is the unit outer surface normal and a. denotes the user supplied surface normal admittance at 
the j-th collocation point. Application of the boundary condition to the incident sound yields the complex source 
vector b g C Nc . 


b j=- 


p 0 c 0 k 


%nc{xj)-nj +a jPinXxj) 


j = l,...,N c 


(15) 


In the FSC v2.0, the overdetermined linear system 

A d = b (16) 

for the unknown monopole source strengths ce C Ns is premultiplied by the complex conjugate matrix (adjoint 
matrix) A^eC^ to produce the N s x N s Hermitian system 

A* Ac = A*b (17) 
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which is then solved using standard Gauss elimination techniques. The preconditioning step, Eq. (17), dominates the 
solution process and requires 0(N c N s ) memory and o(n 2 N s ) multiplications to complete. For large-scale problems 
the computational requirements are well beyond the capabilities of 2008 workstation technology. 

To reduce the computational complexity of solving Eq. (16), sparse, iterative linear algebra techniques based on 
the conjugate gradient method have been incorporated in the FSC v3.0. The methodology features SLFMM 
algorithms for fast matrix- vector multiplication. As a result, the direct matrix-matrix multiply in Eq. (17) is avoided 
and memory requirements are reduced by 95% or more for some of the applications considered here. 

1. Conjugate Gradient Algorithm 

There are several variants of the CGM in the literature (see reference 9 for a review). The algorithm presented 
below was chosen for its compatibility with the FSC/FMM methodology. 

Step 1: For initial guess set r^ = b - Ac^\ p ^ =s ^ = A*r^°\ y 0 =||^| 2 where ||x|| 2 =^|x m | 

m = 1 

Step 2: For A = 0,1,2,... and while > £ for some prescribed tolerance £ compute 


a) 

q^ = 

Ap' z> and a x ■ 

£ 

II 

b) 

c (i+1) 

= c w + a A p w 


c) 


II 

l 

£ 

-b 


d) 


= A V (/L+l) and 

n + i= 

e) 

p (M) 

II 

bi 

+ 

+ 

: 

+ 

^31 



Yx 


Each iteration of the CGM algorithm involves two matrix-vector multiplies. One multiply involves the original 
FSC matrix and the other its adjoint. Direct evaluation of either product requires 0(N c N s ) multiplications and 
constitutes the bulk of the numerical work for the CGM. In the next two subsections, fast matrix- vector multiply 
algorithms are presented that incorporate the SLFMM and reduce the computational complexity an order of 
magnitude. 

Fast Matrix- vector Multiply I 


q = Ap ; 




^Po c o k 


•v+« y 




Sim J : 


■ \ 


(18) 


The sum on the far right in Eq. (18) is exactly of the form of Eq. (2) and therefore can be evaluated by the 
SLFMM in o\ N^\ operations. According to Eq. (12), it is noted that only the portion of the scattering matrix 
involving source-receiver interactions that are in close proximity need to be stored. Additional arrays of size 
o(q 2 N c ) are also required for certain spherical harmonic information that is used repeatedly. 

Fast Matrix- vector Multiply II 

J = A r => ? — r A (19) 


: r 1 


=X/> b 


7=1 


=2/ 

J = 1 


Po c o k 


' V g™ +«/£/ 




( 20 ) 


In Eq. (20), the bar over a variable denotes the complex conjugate. Expanding the gradient operator and 
rearranging terms yields 


3 a N c N c 

=i^-i d '8m + i e j8 l 

1=1 OX l 7=1 7=1 


m = \,...,N s 


( 21 ) 


Where 
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I 


( 22 ) 


d) 


p 0 c 0 k 


fyj 


and 


■ a J r j 


j = 1 = 1,2,3 


Each of the sums over the index j in Eq. (21) is amenable to SLFMM calculations. Again, the computational 
complexity is reduced to j • 

III. Full Scale Applications 

The advantages associated with the computational enhancements implemented into FSC v3.0 are demonstrated 
in this section for three full scale applications. The first two cases highlight usage of the conjugate gradient method 
to obtain the total acoustic pressure at a field of observers bisecting 1) a typical high bypass ratio nacelle, and 2) a 
commercial transport (fuselage, wings, nacelles); the last case features usage of the SLFMM to simulate the 
scattered acoustic pressure field at a plane bisecting the futuristic SAX-40 10 configuration. All calculations presented 
here were performed on a single processor of a Beowulf cluster with sixteen Itanium 2, 1 .4 GHz processors and 32 
GB of shared RAM. 

A. High Bypass Ratio Nacelle 

Because of the large surface area of modern high bypass ratio engine nacelles, predictions using FSC v2.0 are 
limited to frequencies up to 1.5 x BPF using current single-processor computer technology. In order to benchmark 
the CGM implementation with available resources, a nacelle surface similar to the GE-90 was discretized for 1.5 x 
BPF and used with an excitation frequency of 3 x BPF. The incident acoustic field was generated using spinning 
monopoles of unit strength placed at the rotor blade tips. The number of monopoles was chosen so as to excite the 
lowest radiating circumferential mode for this frequency (m = 12). For simplicity, forward Mach number is zero. 
Contours of total acoustic pressure obtained with FSC v2.0 are presented in Figure 2a for a plane bisecting the 
nacelle and Fig. 2b for a plane located 0.05 m forward of the rotor disc. Contours of sound pressure level (SPL) are 
given in Figure 3. Although the computational mesh contains only a fraction of the collocation points (25%) and 
equivalent sources (25%) necessary for proper wave resolution at 3 x BPF, note from the figures that the correct 
trend is simulated - propagation of multiple radial modes through the nacelle openings. 

CGM validation with a null initial guess against known FSC v2.0 results indicated that the method is robust and 
converges rapidly to the correct solution for single component configurations at moderate frequencies. At higher 
frequencies, convergence is improved using a non-zero initial guess obtained with a subset of the required 
equivalent sources. A coarser mesh with point density % of fine was used for this exercise. Instantaneous total 
acoustic pressure contours for the initial guess and final solution are given in Figures 4a and 4b, respectively. Note 
that the CGM fine grid solution is virtually identical to that obtained with the direct method, shown in Fig. 2a. 
Statistics for the different aspects of this application are presented in Table 1. Note that computational time was 
reduced by a factor of 2.5. 



(a) Bisecting plane 


(b) Plane 0.05 m forward of rotor disc 


Figure 2. Instantaneous Re (p’) contours, in Pa, for a full scale GE-90 like nacelle; f = 3 x BPF, 

M = 0.0, m = 12. FSC v2.0 results. 
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(a) Bisecting plane 


(b) Plane 0.05 m forward of rotor disc 


Figure 3. Instantaneous SPL contours, in dB, for a full scale GE-90 like nacelle; f = 3 x BPF, M = 0.0, m = 12. 

FSC v2.0 results. 



(a) FSC-CGM, initial guess 


(b) FSC-CGM, fine grid 


Figure 4. Instantaneous Re (p’) contours, in Pa, for a full scale GE-90 like nacelle; f = 3 x BPF, M = 0.0, 

m = 12. FSC-CGM results. 


Table 1. Statistics for full scale nacelle, f = 3 x BPF (with 1.5 x BPF grid). 


Parameter 

Initial guess 

Fine grid 

Cummulative 

Improvement 
over DM 

No. iterations 

74 

1198 



Solution time, s 

137 

9888 

10052 

2.53 

Direct method (DM) 
solution time, s 

25413 


B. Commercial Transport 

The most challenging application for the FSC with the SLFMM/CGM methodology is high frequency 
predictions for large scale commercial transports. Current FSC practice is limited to excitation frequencies of about 
0.4 x BPF of engine fan noise. In order to test the CGM implementation on an arbitrary, multi-component geometry 
with available computer resources, a configuration similar to the Boeing 111 (wings, fuselage and nacelles) 
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discretized for 0.4 x BPF was used with an excitation frequency of 1 x BPF. Forward velocity was zero. The 
incident acoustic field was generated using the engine noise feedback option, which allows the user to input a pre- 
existing nacelle-only solution into a current FSC run. A nacelle similar to the GE-90 was used for this exercise. 
Because all axial modes are cut-off for 1 x BPF, stationary monopoles at the nacelle centers were used to generate 
the input solution. For simplicity, forward velocity was zero. 

Contours of total acoustic pressure obtained with FSC v2.0 are presented in Figure 5 for a plane bisecting the 
configuration; SPL contours are given in Figure 6. Shielding of engine noise by the wings and fuselage is apparent, 
although the computational mesh contains only a fraction (about 20%) of the required collocation points and 
equivalent sources associated with the frequency under consideration. Results using the CGM algorithm with a non- 
zero initial guess obtained with a % density mesh are presented in Figure 7. Note that the fine mesh solution is very 
similar to that obtained with the direct method (see Fig. 5), even though the iteration limit was reached. Statistics for 
the different aspects of this application are presented in Table 2. For this case, computational time was reduced by a 
factor of 1.38. 


RP 



0.01 

0.008 

0.006 

G.GQ4 

0.002 

0 

- 0.002 

- 0 . 00-1 

- 0.006 

- 0.008 

- 0.01 


Figure 5. Instantaneous contours of Re (p’) for a full scale commercial transport; f = 1 x BPF (with 0.4 x 
BPF grid), M = 0.0. FSC v2.0 results with engine noise feedback. 



Figure 6. Instantaneous SPL contours for a full scale commercial transport; f = 1 x BPF (with 0.4 x BPF 
grid), M = 0.0. FSC v2.0 results with engine noise feedback. 



RP 



(a) FSC-CGM, initial guess 
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(b) FSC-CGM, fine grid 


Figure 7. Instantaneous Re (p’) contours for a full scale commercial transport; f = 1 x BPF (with 0.4 x BPF 
grid), M = 0.0. FSC - CGM results with engine noise feedback. 
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Table 2. Statistics for full scale commercial transport, f = 1 x BPF (with 0.4 x BPF grid). 


Parameter 

Initial guess 

Fine grid 

Cummulative 

Improvement 
over DM 

No. iterations 

466 

2500 a 



Solution time, s 

1090 

26188 

27370 

1.38 

Direct method (DM) 
solution time, s 

37887 


a. Iteration limit reached. 


C. Blended Wing Body Concept (SAX-40) 

The configuration considered in this exercise resulted from the Silent Aircraft Initiative, a collaborative effort 
among industry, academia and government aimed at determining the aerodynamic and performance characteristics 
of a conceptual vehicle using noise reduction as a prime design variable. The final design, known as SAX-40, 
incorporates three boundary layer ingesting engine clusters. A variation of the design, with three pod-type nacelles 
placed above the fuselage, was used here. 

The incident acoustic field, generated using the engine noise feedback option, is that due to spinning monopoles 
of circumferential mode order zero (m = 0) located within the nacelles at the fan face. Available resources limited 
the FSC v2.0 simulations to an excitation frequency of 400 Hz using a subset computational mesh corresponding to 
200 Hz. Contours of SPL (including contributions from all three engines) at a plane bisecting the configuration are 
presented in Figure 8. The plane contains 255 x 510 points. Note from the figure that, although the acoustic waves 
are not properly resolved for the frequency of interest, shielding of fan inlet noise by the body is apparent. Contours 
of scattered acoustic pressure at the same observer field are given in Figure 9a. The SLFMM implementation, in 
conjunction with equivalent source strengths previously calculated with FSC v2.0, was used to simulate the scattered 
acoustic field presented in Figure 9b. Observe that the fields are virtually identical. Statistics for the different aspects 
of this test are included in Table 3. Computational time was reduced by a factor of 3.78. 




SPL: 60 62 64 66 68 70 72 74 76 78 80 82 84 86 88 90 


Figure 9. Instantaneous SPL contours, in dB, for the SAX-40 conceptual design with pod nacelles; 
f = 400 Hz (with 200 Hz grid), M = 0.0, m = 0. FSC v2.0 results. 
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(a) FSC v2.0 (b) FSC - SLFMM 

Figure 10. Instantaneous Re (p ' scat ) contours, in Pa, for the SAX-40 design with pod nacelles; f = 400 Hz 

(with 200 Hz grid), M = 0.0, m = 0. 


Table 3. Statistics for SAX-40 configuration, f = 400 Hz (with 200 Hz grid). 


Parameter 

Value 

Improvement over 
DM 

Terms in spherical expansion 

23 


Non-empty boxes 

325 


Average points per box 

400 


Box size, m 

2.165 


Storage (relative to DM), % 

3.5 


Time, s 

2399 

3.78 

Direct method time, s 

9072 


IV. Concluding Remarks 

The NASA Fast Scattering Code is a versatile analysis tool developed to conduct aeroacoustic noise reduction 
studies. The program employs the equivalent source method to solve a three-dimensional exterior Helmholtz 
boundary value problem. FSC v2.0 uses dense matrix solution techniques that are memory intensive and time 
consuming, limiting applications to problems of moderate size. Calculations are dominated by the matrix-matrix 
multiplication required by the pre-conditioning step of the solution process. Application of the single-level fast 
multipole method and sparse linear algebra techniques eliminates the matrix-matrix multiply and has led to an order- 
of-magnitude reduction in computational resource requirements, thereby doubling the excitation frequency 
achievable for a given configuration. In the program FSC v3.0, iterative conjugate gradient methods have been 
coupled to fast matrix-vector multiplication algorithms based on the SLFMM. At the end of the solution process, the 
SLFMM is deployed again for the rapid summation of the acoustic field. Sample calculations for several 
configurations of importance have been presented to demonstrate the SLFMM analysis and speed-up gains. 

Further computational advances are possible with the SLFMM technology as presented here, and will be 
addressed by the authors in future releases of the FSC. First, o(q 3 ) translation algorithms exist in the literature for 
faster evaluation of the regular expansion coefficients in Eq. (11) which is an o[q A ) algorithm. Second, multi-level 
fast multipole methods (MLFMM) offer an additional degree of computational time reduction. By repeatedly 
subdividing the computational boxes into smaller boxes, analysis shows that the complexity of the SLFMM can be 
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reduced from o\ N^ 3 ] to o(viogV) by application of the MLFMM algorithms 8 . These two research aspects 
account for marginal ^gains in the FSC performance. More profound computational advantages can be obtained 
through the use of multi-processor computer technology. Step 2 in the evaluation of the SLFMM algorithm involves 
the calculation of the regular expansion for all points in each receiver box. Since receiver box calculations are 
independent of each other, it is possible to distribute the computational work for the boxes over an array of 
processors reducing the work load by a multiplicative factor approaching the number of processors. Incorporation of 
these additional research methods is expected to be completed by early 2009. 
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